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ABSTRACT 

We propose a model for stellar binary systems consisting of a magnetic and a 
non-magnetic white-dwarf pair which is powered principally by electrical energy. In 
our model the luminosity is caused by resistive heating of the stellar atmospheres 
due to induced currents driven within the binary. This process is reminiscent of the 
Jupiter-Io system, but greatly increased in power because of the larger companion 
and stronger magnetic field of the primary. Electrical power is an alternative stellar 
luminosity source, following on from nuclear fusion and accretion. We find that this 
source of heating is sufficient to account for the observed X-ray luminosity of the 9.5- 
min binary RX J1914+24, and provides an explanation for its puzzling characteristics. 
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1 INTRODUCTION 

It has been observed directly by the Hubble Space Telescope 
that the movement of Io through Jupiter's magnetic field 
causes heating in the Jovian atmosphere (Clarke et al. 1996) . 
This is because a conducting body transversing a magnetic 
field produces an induced electric field. When the circuit is 
closed, a current will be set up, resulting in resistive dissipa- 
tion. The Jupiter-Io system therefore operates as a unipolar 
inductor (Paddington & Drake 1968; Goldreich & Lynden- 
Bell 1969) . Another potential cosmic unipolar inductor could 
be a planet orbiting around a magnetic white dwarf (Li, Fer- 
rario & Wickramasinghe 1998). These systems have a similar 
configuration, with the differences being their orbital period 
and separation, the masses and radii of the two components, 
and the magnetic moment of the magnetic object. 

We propose that binary stars consisting of a magnetic 
and a non-magnetic white dwarf can also be cosmic unipo- 
lar inductors (Fig. 1). Close binaries of this type can have 
short periods and secondaries larger than planet-sized bod- 
ies. Provided that the spin of the magnetic white dwarf and 
the orbital rotation are not synchronised (so that the sec- 
ondary is in motion relative to the magnetic field) and that 
the density of the plasma between the white dwarfs is high 
enough, unipolar induction will operate efficiently. 

Gravitational waves carry away the orbital angular mo- 
mentum efficiently from a short-period (<2 hr) binary sys- 
tem but not the stellar spin momenta directly. In a mag- 
netic and non-magnetic white-dwarf pair with only the non- 
magnetic star tidally locked, the magnetic star will be spun 



up retrogradely in the orbital rest frame as the binary or- 
bit shrinks. Alternatively, the magnetic white dwarf could 
be spun up by accretion in a previous epoch in which mass 
transfer occurred, so that it has a spin faster than the orbital 
rotation. Asynchronous rotation can therefore occur as long 
as a stellar component is not tidally or otherwise locked. 

The dissipative power of a white dwarf-planet pair with 
an orbital period P — 10 hr is estimated to be ~ 10 29 erg s _1 
(Li, Ferrario & Wickramasinghe 1998). The output of a bi- 
nary white-dwarf pair with a period of tens of minutes will 
be significantly higher and should have detectable observa- 
tional consequences. 

In this paper, a simple model for a unipolar inductor 
consisting of a magnetic and a non-magnetic white-dwarf 
pair in a close orbit is presented. The basic features of 
these binaries and their observational properties are dis- 
cussed. We propose that the short-period soft X-ray source 
RX J1914+24 (Cropper et al. 1998; Ramsay et al. 2000a) is 
a candidate unipolar inductor consisting of a magnetic and 
non-magnetic white-dwarf pair. 



2 WHITE-DWARF PAIRS AS COSMIC 
UNIPOLAR INDUCTORS 

When a non-magnetic conductor of linear size R transvers- 
ing a magnetic field B with a velocity v, the induced e.m.f. 
across the conductor is "I> ~ R\E\, where E = (v x B)/c (c 
is the speed of light). Thus, the corresponding e.m.f. across 
a non-magnetic white dwarf in orbit with a magnetic white 
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Figure 1. A schematic illustration of a unipolar inductor consist- 
ing of a magnetic and non-magnetic white-dwarf pair in a close 
binary orbit. 
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where G is the gravitational constant, q (= M2/M1) is the 
ratio of the non-magnetic to the magnetic white dwarf mass, 
R2 and Ri are the radii of the non-magnetic and the mag- 
netic white dwarf respectively, /ii is the magnetic moment of 
the magnetic white dwarf. The degree of asynchronism a is 
defined as the ratio of the spin angular speed of the magnetic 
white dwarf ui to the orbital angular speed lu (= 2n/P). 
(Here and elsewhere, we consider that the ant i- clockwise di- 
rection is positive.) We show in Figure 2 the induced e.m.f. 
for different system parameters. The induced e.m.f. depends 
strongly on the orbital period, the degree of spin-orbit syn- 
chronism, and the mass (radius) of the non-magnetic white 
dwarf. The dependence of the mass of the magnetic white 
dwarf is weaker. 

If the space between the white dwarfs is filled with 
plasma, the induced e.m.f. will drive currents along the 
magnetic field lines connecting the two white dwarfs. Be- 
cause there is substantial resistance in the white-dwarf at- 
mosphere, the white dwarfs act as the dissipative compo- 
nents in this electrical circuit. The total electrical power 
dissipation in the two white dwarfs is 



W = I 2 (TZ 1 +TZ 2 ) 
= $ 2 /(^i+^i), 



(2) 



where I is the total current, and 1Zi and TI2 are the effec- 
tive resistance of the magnetic and the non-magnetic white 
dwarfs respectively. An object with a length L and an area 
A has a resistance 1Z = L/Aa (where a is conductivity). 
Therefore, the ratio of the effective resistances of the white 
dwarfs is 
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Figure 2. The induced e.m.f. across the non-magnetic secondary 
white dwarf as a function of orbital period for binaries with vari- 
ous parameters. Solid lines a, b and c correspond to systems with 
a primary magnetic white dwarf of 1.0 Mq with a spin 1 part 
of 1000 deviated from synchronous rotation with the orbit. The 
masses of the secondaries in these systems are 0.1, 0,5 and 1.0 Mq 
respectively. The dotted line corresponds to systems consisting 
of a 0.1-Mq non-magnetic white dwarf, and a 0.7-Mq magnetic 
white dwarf with 1 part of 1000 off spin-orbit synchronism; and 
the dashed line corresponds to similar systems but with a 1.3- 
Mq primary magnetic white dwarf. Systems with magnetic and 
non-magnetic white-dwarf pairs of mass 1.0 and 0.5 Mq and spin- 
orbital asynchronism of 1 part in 100 for the magnetic white dwarf 
is represented by the dot-dashed line. The magnetic moments of 
the systems are fixed to be 10 a2 G cm 3 . 



where <ti and 02 are the conductivity of the magnetic and 
the non-magnetic white dwarf respectively, Ah\_ and A/12 are 
the thickness of the dissipative surface layers of the white 
dwarfs, and / is the fractional effective area of the magnetic 
poles (hot spots) on the surface of the magnetic white dwarf. 
As / < 1 (see Appendix A), the effective resistance of the 
magnetic white dwarf is significantly larger than the non- 
magnetic white dwarf. 

As the currents pass through both white dwarfs, the 
ratio of the power dissipation in the magnetic primary to 
that of the non-magnetic secondary is W1/W2 = TZi/Ti.2. 
Taking account of the geometry of the current loops, we 
obtain 
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(see Appendices A and B for details), where AR 2 is the 
thickness of the secondary's atmosphere and (5 is a structure 
factor of the order of unity. The factor J(e) depends of the 
radii of the white dwarfs relative to the orbital separation. 
Its value is of the order of unity for white-dwarf pairs with 
orbital periods less than an hour. 

The conductivity of plasma of an electron temperature 
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T e is given by 
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(Spitzer & Harm 1953) where k is the Boltzmann constant, 
T c is the electron temperature, m c is the electron mass, e 
is the electron charge, Z is the ion charge number, and In A 
is the Coulomb logarithm. The factor 7 depends on the ion 
charge number Z, which has values between 0.6 (Z = 1) and 
1 (Z — > 00) (see Alfven & Falthammar 1963). For a white- 
dwarf atmosphere with T c ~ 10 5 K, the conductivity a ~ 
10 13 - 10 14 esu. Since the conductivities of the atmospheres 
of white dwarfs are similar to each other, the majority of 
the electrical power will be dissipated in small regions at the 
footpoints of the current-carrying field lines on the surface 
of the magnetic white dwarf. 



3 ENERGY AND ANGULAR-MOMENTUM 
CONSERVATION 

3.1 Power dissipation 

If the degree of spin-orbit asynchronism (1 — a) is specified, 
the total power dissipation of the binary can be calculated 
directly using equations (1), (2), (4) and (5). We show in 
Figure 3 the total power dissipation W as a function of the 
orbital period P, for systems with different primary and sec- 
ondary white-dwarf masses Mi and M2. The general trend 
of the dependence of W on the system parameters is as 
follows: (i) the power dissipation W is larger for larger de- 
gree of spin-orbit asynchronism; (ii) W increases when P 
decreases; (iii) the larger is the radius of the non-magnetic 
white dwarf R2, the larger W will be; and (iv) W decreases 
slightly when M\ increases. Moreover, provided that there 
is a slight spin-orbit asynchronism (1 part in 1000 is suffi- 
cient) and that the asynchronously rotating white dwarf has 
a moderate magnetic moment (~ 10 32 G cm 3 ), the electrical 
power dissipation of the systems can reach the solar luminos- 
ity Lq (3.9 x 10 33 erg s _1 ), when the orbital period is short 
enough (S;10 min). Clearly, the electrical power generated by 
these close binaries is substantially larger than the expected 
intrinsic luminosities of the white-dwarf components, which 
is <C Lq . 

It is worth noting that most of the luminosity origi- 
nates from the two hot spots near the magnetic poles at 
the surface of the magnetic white-dwarf. If the footpoints 
of the current-carrying fieldlines remain fixed with respect 
to the secondary, as in the Jupiter-Io system (Clarke et 
al. 1996), asynchronous rotation will not manifest itself in 
orbital-period changes, but will possibly be visible through 
changes in the spot latitudes if the dipole is not aligned with 
the spin axis. 



3.2 Spin-orbit coupling 

The consequence of the energy dissipation is the creation of 
a Lorentz torque, which tends to synchronise the spin of the 
magnetic white dwarf with the spin of the non-magnetic 
white dwarf, and hence with the orbital rotation. Thus, 
through spin-orbit coupling, energy and angular momentum 
can also be exchanged between orbit and the white dwarfs. 




Figure 3. The total electrical power dissipation as a function of 
orbital period. The degree of deviation from spin-orbit asynchro- 
nism (I — a) is 1 part in 1000 in the top panel and 1 part in 100 in 
the bottom panel. As in Fig. 2 the solid lines a, b and c correspond 
to systems with a non-magnetic white-dwarf secondary of 0.1, 0.5 
and 1.0 Mq. The mass of the magnetic primary is 1.0 Mq. The 
dotted line corresponds to systems consisting of a 0.7-Mq mag- 
netic white dwarf and a 0.1-Mq non- magnetic white dwarf, and 
the dashed line corresponds to systems consisting of a 1.3-Mq 
magnetic white dwarf and a 0.1-Mq non-magnetic white dwarf. 
The magnetic moment of the primary is fixed to be 10 32 G cm 3 . 



If the secondary white dwarf is tidally locked to the orbit, as 
assumed in our model, the evolution of the spin of the mag- 
netic primary white dwarf wi, the orbital rotation uj , and 
the degree of synchronisation is governed by the equations: 
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(see Appendices C, D and E), where 
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are expected to be rare in comparison with 'ordinary' white- 
dwarf pairs. 




Figure 4. The total electrical power dissipation (solid lines) and 
the power of gravitation radiation emitted from the binary (dot- 
ted lines) as a function of mass of the non-magnetic white dwarf. 
The orbital period is fixed to be 9.5 min. Curves a, b, c and 
d correspond to systems with a 0.5, 0.7, 1.0, and 1.3-Mq mag- 
netic white dwarf. The magnetic moments of the white dwarfs 
are 10 32 G cm 3 . The spin of the magnetic white dwarf is 1 part 
in 1000 deviated from synchronous rotation with the orbit, i.e. 
(1 - a ) = 0.001. 



+ «)/(<*>) 
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The term E cyi is the energy loss due to external process 
(such as gravitational radiation). The factor f(io ) can be 
expressed in terms of orbital parameters, which is simply 
(R 2 /a) 2 . Clearly, the spin-orbit evolution is jointly deter- 
mined by the rate of energy loss due to gravitation radiation 
E gI and the electrical power dissipation W. 

For short-period binaries gravitational radiation is an 
efficient process to extract energy from the binary orbit 
(Paczhyski 1967; Faulkner 1971). As an illustration we show 
in Figure 4 the electric power dissipation W and the rate 
of energy loss from the orbit due to gravitation radiation 
E ST for binaries with an orbital period < 10 min and spin- 
orbital asynchronism of 1 part in 1000. The power loss due to 
gravitational radiation is generally larger than the electrical 
power dissipation, except for systems with a very low-mass 
non-magnetic white dwarf. 

As W oc (1 — a) 2 > (from Equations 1 and 2), energy 
is extracted from the spin of the magnetic white dwarf, if it 
is a fast spinner a > 1 (Equation 9). When a < 1, energy is 
injected to spin up the magnetic white dwarf. Because both 
E gI and g(i0 o ) are negative, when energy is extracted from 
the orbit to spin up the magnetic white dwarf, the orbit 
shrinks and the orbital angular speed increases. For white- 
dwarf pairs consisting of a l.O-M© magnetic primary and a 
0.1-Mq secondary with a 10-min orbit and a spin-orbital 
asynchronism of 1 part in 1000, the typical evolutionary 
timescale for a is about 10 3 years. As the synchronisation 
timescale due to unipolar induction is significantly shorter 
than the timescale of the gravitational-radiation power loss, 
which is about 10 6 years for these systems, these systems 



3.3 Effects of the induced magnetic field 

In the unipolar-inductor model that we consider above, the 
magnetic field generated by the currents between the white 
dwarfs has been neglected. Because of this induced magnetic 
field, the magnetic-field configuration is no longer dipolar as 
for the assumed intrinsic field of the primary white dwarf. 
As a result, the currents and the magnetic field lines are not 
aligned, and there will an additional drag to the motion of 
the secondary white dwarf. This plasma-inertia effect can be 
non-negligible. 

The importance of plasma-inertia effects in the Jupiter- 
Io system was investigated by Drell, Foley & Ruderman 
(1965) in the linear approximation and by Neubauer (1980) 
using a more appropriate non-linear treatment. These stud- 
ies showed that the non-alignment of the currents and mag- 
netic field is associated with the Alfven waves — called 
Alfven wings — standing in the frame of Io, causing dis- 
sipation in the polar regions of Jupiter as well as in Io. 

Here in this white-dwarf pair, the standing Alfven waves 
are in the frame of the non-magnetic secondary white dwarf. 
The significance of the plasma-inertia effects can be ex- 
pressed by the Alfven Mach number 



M A = — , 

va 



(12) 



where 5v is velocity of secondary white dwarf related to the 
magnetosphere of the primary magnetic white dwarf, given 
by 5v ~ (1 — a)-iva/P w 10 9 (1 — a) cm s _1 , and va is the 
Alfven speed: 



v A = T5 x 10 J 
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(n e is the electron number density in the plasma). For the 
parameters of our interest, Ma 1 is generally satisfied, 
and the plasma-inertia effects can be considered as a per- 
turbation. The results that we have obtained are therefore 
generally valid. 

The large currents in the system can also cause signifi- 
cant distortion of the magnetic field near the primary mag- 
netic white dwarf. The geometry of the energy-dissipation 
region and the electrical circuit of the system are therefore 
not as simple as described in our idealised model. A particu- 
lar issue is whether or not steady currents can be maintained 
given the presence of a large j x B force, when the current 
and the magnetic field are not aligned. This could eventually 
lead to the breakdown of the circuit and quench the opera- 
tion of the unipolar induction. Nevertheless, if the unipolar 
induction does operate, in spite of the non-alignment of cur- 
rent and field, the energy and angular-momentum budget 
should be similar to that of the unipolar-inductor model 
that we have presented above. (Detail investigation of the 
effects due to severe field distortion is beyond the scope of 
this paper, and will be left for future study.) 
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4 RX J1914+24: A POSSIBLE UNIPOLAR 
INDUCTOR? 

We propose that the short-period binary system RX 
J1914+24 is electrically powered. This system emits soft X- 
rays which are modulated on a period of 9.5 min, and whose 
folded light curve is that expected from emission originating 
from one or two small spots. The optical/infra-red flux is 
also modulated only on the 9.5-min period but maximum 
light occurs ~ 0.4 orbital cycles earlier than the soft X-rays 
(Ramsay et al. 2000b). 

It has been argued that the system is a magnetic 
white dwarf-white dwarf binary, which is synchronised by 
the strong magnetic interaction between the white dwarfs 
(Cropper et al. 1998). An accretion scenario has been pro- 
posed (Ramsay et al. 2000a), in which the optical flux is 
due to irradiation by the accreting material, and which can 
largely account for the observations. However, to explain 
the lack of optical polarisation (Ramsay et al. 2000a), the 
magnetic-field strength on the primary must either be too 
high or too low to be detected in the optical/infra-red. Also, 
a major concern is that the accretion model is unable to ex- 
plain the observed lack of strong emission lines (Ramsay et 
al. 2000b) which result from heating of the accretion stream 
in mass-transfer systems. 

If the secondary in RX J1914+24 fills its Roche lobe 
(Cropper et al. 1998), M 2 ~ O.1M . A Roche-lobe-filling 
star is more easily tidally locked into synchronous rotation 
with the orbit. In the unipolar-inductor model proposed 
above, the secondary is allowed to lie within its Roche lobe, 
and so this mass is a lower limit. The mass of the primary 
is, however, unconstrained. Based on ROSAT and ASCA 
measurements and an estimated distance of 200-500 pc, the 
deduced luminosity of RX J1914+24 is in the range 4 x 10 33 
to 1 x 10 35 erg s _1 (Ramsay et al. 2000b). Provided that 
the spin-orbit asynchronism of the magnetic white dwarf is 
about 1 part in 1000 or larger, the resistive heating is suffi- 
cient to power these luminosities (See Figures 3 and 4). From 
Equation (4), the resistive dissipation power is estimated to 
be ~ 10 31 — 10 32 erg s _1 in the secondary's atmosphere. 
The result is less sensitive to other system parameters such 
as the mass of the magnetic white dwarf. A comparison be- 
tween the observed and predicted luminosities suggests that 
the secondary is 0.15 M Q , and thus close to filling its Roche 
lobe and tidally locked. 

The predicted area of each hotspot on the primary is 
~8x 10 14 cm 2 (from Appendix A) . If we use the observed 
luminosity of 10 34 erg s _1 (assuming a distance d — 300 pc), 
we obtain a blackbody temperature of 50 eV, consistent with 
the measured value of 55 eV from the ROSAT and ASCA 
spectral fits (Ramsay et al. 2000b). This consistent area is 
corroborating evidence that the underlying geometry of the 
model is valid. 

For a system with a l.O-M© primary, a O.l-M© sec- 
ondary and an orbital period of 9.5 min, less than 1 percent 
of the X-ray flux emitted from the hotspots primary is in- 
tercepted by the secondary. In equilibrium, the intercepted 
flux will be reradiated, resulting in an increase in the lumi- 
nosity of the irradiated hemisphere of the secondary. The 
additional luminosity arising from irradiation heating is es- 
timated to be ;$10 32 erg s _1 . Thus, the power of irradiative 
heating of the secondary is similar to the power of electrical 



heating. Given the fact that the two heating processes have 
similar efficiency, one would not expect a large temperature 
difference in the irradiated and the unirradiated hemisphere 
of the secondary. This is consistent with the observation 
that the amplitude of variation in the optical luminosity in 
RX J1914+24 is small (Ramsay et al. 2000a). Moreover, as 
the power of irradiative heating is of the same order of the 
energy flux from below the atmosphere due to resistive heat- 
ing, the presence of moderate heat conduction is sufficient 
to maintain an approximate isothermal layer down to unit 
optical depth. Because of the absence of an temperature in- 
version layer, the secondary will not show emission lines in 
its spectrum: this is consistent with the lack of prominent 
emission lines in the optical spectrum (Ramsay et al. 2000b) . 



5 SUMMARY 

We propose that short-period magnetic and non-magnetic 
white-dwarf pairs with short orbital periods (~ 10 min) are 
efficient cosmic unipolar inductors. Provided that the spin 
of the magnetic component and the orbit are not in perfect 
synchronism, a large e.m.f. can be produced across the non- 
magnetic white dwarf. The resistive dissipation in the white 
dwarfs is sufficient to power luminosities significantly above 
solar values; most power is dissipated at the hot spots on the 
surface of the magnetic white dwarfs, which are footpoints 
of the field lines connecting the two stars. Electrical power 
is therefore an alternative luminosity source, following on 
from nuclear fusion and accretion. 

The X-ray source RX J1914+24 is a candidate unipolar 
inductors consisting of a magnetic and non-magnetic white- 
dwarf pair. The two small X-ray spots on the magnetic white 
dwarf predicted by the unipolar-inductor model are compat- 
ible with the X-ray light curve of RX J1914+24. The lumi- 
nosity and temperature predicted by the model is also in 
agreement with the observed values derived from fits to the 
X-ray spectra. The model also explains the variation in the 
optical/infra-red luminosity, and the detection of only a sin- 
gle period. The variations in the long term X-ray intensity 
can be attributed to variations in the current flow. The two 
main inadequacies of the current accretion model (Ramsay 
ct al. 2000a) — the lack of any polarised flux and the lack 
of any detectable line emission (Ramsay et al. 2000b) arc 
naturally explained. 
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APPENDIX A: 
SPOTS 



GEOMETRY OF THE HOT 



Consider a system consisting of a magnetic and a non- 
magnetic white dwarf in close orbit with a separation a. 
The magnetic white dwarf has a dipolar field, with the fam- 
ily of field lines given by r = C sin 2 #, where C is a constant 
labeling different field lines. Thus, the geometry of the pole 
is determined by the field lines that connected the two white 
dwarfs, and it is specified by the parameter a and 6, where 
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(see Figure 5, top and second panels). 

If the radius of the stars Ri and Ri are significantly 
smaller than the orbital separation a, then a w 7? 2 (i?i/a) 3 / 2 
and b ~ (_R 2 /2)(-Ri/a) 3 ' 2 . It follows that b/a is approxi- 
mately 0.5 for a wide range of orbital parameters. Moreover, 
the fractional area of the dissipative region on the white- 
dwarf surface / < (ab/Rf) <C 1. 



APPENDIX B: RESISTIVE DISSIPATION IN 
THE TWO WHITE DWARFS 

For the non-magnetic white dwarf, the conductivity of the 
core ct co is much larger than the conductivity of the atmo- 
sphere <J2 (i.e. cr co 3> 02), so that most resistivity dissipa- 
tion occurs at the atmospheric layer (See Figure 5, bottom 
panel). The rate of energy dissipation is therefore 
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Figure 5. (Top and second panels) The top and side elevation 
views of the system showing how maximum height and width 
of the hot spot (26 and 2d respectively) at the surface of the 
magnetic white dwarf are determined. (Third penal) Schematic 
illustration of the current sheets in the atmosphere of the mag- 
netic primary white dwarf. At depth H below the white-dwarf, 
collisions are efficient enough to allow the current to cross the 
magnetic-field lines. Bottom panel) Schematic illustration of the 
current flow in the non-magnetic secondary white dwarf. 
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If the current density J is roughly homogeneous, then 
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W 2 = 4ttRIAR 2 — 
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4 I 



(17) 



The exact values of H and Ad are uncertain. Li, Wickra- 
masinghe & Ferrario (1998) argued that roughly H ~ Ad. 
As 02 ~ ffi, nJ(e)/8 ~ 0(1) and (i? 2 / AR 2 )(a/ Ri) :i/2 > 1, 
most of the electric energy will be dissipated at the two hot 
polar spots of the magnetic white dwarf. 



where I — n R2J2 is the current flowing across the star. 

For the magnetic white dwarf, the dissipation at each 
polar spot is 



J spot 



(18) 



It is not trivial to determine the functional form of the 
current density Ji(x) at the pole. For simplicity, we follow 
Li, Wickramasinghe and Ferrario (1998) and consider an 
approximation that the current density is non-zero in two 
arc-like regions. The thickness of the arc is Ad, and the 
length and width of the region are 25 = Ri sin 9A<j> and 
2b = R1A8, where 8 is the colatitude of the centre of the hot 
spot. The current density is given by Ji = Ji sin£ (Figure 
5, third panel). We assume that at a depth ff(<C Ri) the 
Pendersen conductivity (see Alfven & Falthammar 1963) 
becomes sufficient large to allow the current to cross the 
field lines and close the circuits. 

The current flowing through the spot is therefore 



dxdy Ji 



/ ■ 

J spot 

/■tt/2 

Jio(4fiAd) / 1 
Jo 



rfC Qi(£;e) 



and the resistive dissipation is 

j2 fir/2 

w (i) = ^L^aHAd) / 
°"i Jo 
The function Q n is given by 

Qn(£; e ) = sin n ^-\/l - e 2 cos 2 £ 



(19) 



(20) 



(21) 



where e = yl — (b/a) 2 . 

Using the expression for a in Appendix A, we obtain 
the rate of total dissipation at the magnetic white dwarf: 

2W^ 



Wi 



IP 

2 CTi 



H 
Ad 



3/2 



R 2 



(22) 



where 



J(e) = 



[j; /2 dtQ^e)\- 



E 



(2n- 1)!! 



(2n)! 



e z H — sm 

e 



(n + l)(2n-l) 
-2 



(23) 



The ratio of the dissipation rates at the two white dwarfs is 
therefore 



Wi 
W 2 



;J(e) 



H 
Ad 



R2 
AR 2 



3/2 



(24) 



APPENDIX C: ENERGY CONSERVATION 

Consider a system with two gravitationally bound, rotating 
objects, with moments of inertia I\ and 7 2 , in circular motion 
around each other. Let the angular speeds of their rotation 
be oj\ and UJ2 respectively. The total mechanical energy of 
the system is 



E = ±I lLU 2 1 + ll2co 2 + (T+V) 



(25) 



T and V are the kinetic and potential energy of the orbit, 
which are related by 



T = - I u 

1 GMlq 
~ 2 a 



(26) 



where the effective moment of inertia of the orbit I — 
Mia 2 \q/{1 + q)] and the orbital angular speed uj — 
[GM,(l + ? )/a 3 ] 1/2 . 

If the spin of the second object is in perfect synchronism 
with the orbit (i.e. L02 = u> ), then 



E 



-hui 1 



Mi 



1 + 9 



':M iq R 2 2 



1 T 2 1 GMlq 
2 hUJl ~2^- 



a - -{l + q)Rl 




(27) 



(Here, we have assumed that the objects are spherical, so 
that their moments of inertia is 2/5 of their mass times the 
square of their radius.) It follows that 



E = Jiwiwi - \;GMlq^- 
2 dt 



1 2 /-, -,Rl 
7-7(l + 9)-T 



= I 1 LU 1 UJ 1 + 
= IlUlldll — 



1 GMlq 



r 2 



:(l + q)f(w ) 



2 Co 

3 U! 



,(28) 



where /(w ) = {i?|oJo/[GMi(l+g)]} 2/3 . If the second object 
is close to filling its Roche-lobe, then R2 » oA[g/(l + 9)] 1 '' 3 , 
and 



1 2 

E JlWlWl + 2^ o0; o 



;A^(l+ 9 )^9 



1/3 2/3 



(29) 



where A = 0.462. Clearly, when 9 < 0.7, A 2 (l + g) 1/3 g 2/3 < 
5/6, and hence the square bracketed term is positive. 
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APPENDIX D: WORK DONE BY THE 
DISSIPATIVE TORQUE 

The torque that accelerates/decelerates the spin of the ob- 
ject 1 is 

r = (30) 

and the torque that changes the spin of object 2 and the 
orbital rotation is 

- t d - r d, _ _ d . . 

r cx - T = 1 — UJ + I2-7-LO2 + LOo-7- lo , (31) 

at at at 

where r cx is an external torque. 

From the definition of the synchronisation a = lui/lu , 
we have 



a 
a 



Wl 



L0 o 
LJ 



(32) 



The perfect-synchronism condition for object 2 and the orbit 
implies that 



r cx - r = (7 + 7 2 )w + u Io ■ 



(33) 



Thus, the timescale for the change in the orbital angular 
speed is 



U - = K 

0J o 



ali 



lo 



(Jo + I 2 )lu \ lo + h J a {Io + h) 
and the timescale of spin evolution of object 1 is 

Tex . a Ic 



^ = K 

LOl 



a 

+ - 



(lo + h)u)o a (lo + h) 



, (34) 



(35) 



where K = [l + a/i/(/ + /2)]~ 1 . The first term in the square 
of above expression is due to the synchronisation process, the 
second term is the contribution of the external torque, and 
the last term is caused by the readjustment of the orbital 
separation. 

The work done by the (dissipative) synchronisation 
torque is 



W = l-F- (oJ — t3i)| 
= |/iwiw (l - a) I 

7iw„a(l - a) 



a | Up 
a u) 



(36) 



It is straightforward to show that the work done can also be 
expressed as follows: 



where 
fl(w ) = 



(1 + 9) 



G 2 M\u 2 



1/3 



1- g(l + 9)/(w ) 



(39) 



The functional form f(ui ) can be found in Appendix C. For 
dissipative synchronisation torque proportional to (lu — cUi), 
we have 



E di! 



-W 

-(1 — u)I 1 uj u 1 



(40) 



It can be shown that the evolution of the spin of the mag- 
netic white dwarf and the orbital rotation are determined 
by the following equations: 



L0 o 

a 
a 



1 



w 



(1-a) 
W 



(1-a) 



1 + 



ahojo 



(41) 
(42) 



W = 



hu) a(l - a) 



1 + 



ah 



a 

+ - 



lo 



(lo + h)u> a (J + h) 



(37) 



APPENDIX E: SPIN-ORBIT EVOLUTION 

Suppose the energy loss is represented by a dissipative term 
due to synchronisation torque E A - lss and an additional term 
TJex (e.g. due to gravitational radiation) . Then we have 

-Ediss — E cx = /iwiwi + g{u) ) ( "J" J ' ( 38 ) 



© 0000 RAS, MNRAS 000, 000-000 



